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Abstract: We investigate the collider phenomenology of gluino-bino co-annihilation scenario 
both at the Tevatron and 7 TeV LHC. This scenario can be realized, for example, in a class 
of realistic supersymmetric models with non-universal gaugino masses and t — b — r Yukawa 
unification. The NLSP gluino and LSP bino should be nearly degenerate in mass, so that the 
typical gluino search channels involving leptons or hard jets are not available. Consequently, 
the gluino can be lighter than various bounds on its mass from direct searches. We propose 
a new search for NLSP gluino involving multi-b final states, arising from the three-body 
decay g — > bt>Xi- We identify two realistic models with gluino mass of around 300 GeV for 
which the three-body decay is dominant, and show that a 4.5 a observation sensitivity can be 
achieved at the Tevatron with an integrated luminosity of 10 fb _1 . For the 7 TeV LHC with 
50 pb _1 of integrated luminosity, the number of signal events for the two models is O(10), to 
be compared with negligible SM background events. 
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1. Introduction 

Low-scale super symmetry (SUSY) is arguably the leading candidate of new physics and will 
be seriously tested at the CERN Large Hadron Collider (LHC). It provides an elegant solution 
to the gauge hierarchy problem and with conserved imparity, the thermal relic abundance of 
the lightest neutralino (LSP) can often be identified with dark matter, consistent with the 
current cosmological observations |j|] . The recent results reported by the CDMS-II experiment 
may indicate the presence of dark matter with mass of around 0(100) GeV Q. To reproduce 
the required thermal relic abundance, a pure wino or Higgsino dark matter should have mass 
of a few TeV due to rapid annihilation rate. On the contrary, the annihilation rate of a pure 
bino with mass of around 100 GeV is too slow, and leads to excessive relic abundance ||. 
A variety of constraints from low energy flavor physics and CP violating physics typically 
favor scenarios with heavy sfermions [^, ||, ||. Consequently, pure bino self annihilation is 
further suppressed if scalars in the t-channel become heavier. To enhance the annihilation 
rate of pure bino, there are generally two categories of models available, and both will lead 
to interesting phenomenology at the LHC. In one scenario the dark matter is a bino- wino or 
bino-Higgsino mixture Q, so that XiXi annihilation is enhanced via the enlarged coupling. 
The mass degeneracy Mi ~ M2 leads to a nearly degenerate chargino-neutralino spectrum ||. 
The second scenario is co-annihilation || [l^] , and if the scalars are heavy, the required relic 
abundance can be achieved through co-annihilation with a strongly-interacting particle such 
as the gluino, namely Xi9,99 ~^ ff @> 0> 12 > [L§- This scenario can be realized in a realistic 



class of t — b — t Yukawa unified models [11, [L2, 14]. In the gluino-bino co-annihilation 



scenarios, the gluino is the next-to-lightest-supersymmetric-particle (NLSP), and the mass 
splitting betwen g and Xi 1S relatively small, namely 0] 



Note that the Sommerfeld enhancement of bino-gluino co-annihilation cross section is rather 
mild, such that Eq. ((□]) is altered by only 2 - 3% p|, p . 

The relations between gaugino masses will be crucial in understanding the nature of su- 
persymmetry breaking and of the underlying theory at ultra-high energy scale. For instance, 
supersymmetric grand unification, string dilaton SUSY breaking, and minimal gauge medi- 
ation all predict the gaugino mass relations M 3 /g% = M 2 /#f = M x jg\ . With NLSP gluino, 
one must invoke non- universal gaugino masses at Mqut- This not only implies very differ- 
ent physics from the above models but can also give rise to very different phenomenology at 
hadron colliders. We focus in this paper on the phenomenological implications of bino-gluino 
co-annihilation scenarios at hadron colliders. 

Being a color octet fermion, gluino pair production is the most promising discovery 
channel for super symmetry at hadron colliders. In an environment such as this with huge 
QCD jet backgrounds, isolated charged leptons ((J^, e^) and 6-jets usually play an important 
role in the searches. For the most well studied scenarios where the charginos are lighter than 
gluinos, the Majorana nature of the gluino will result in same-sign chargino signature (jets 
plus X^X* 1 )- This eventually leads to same-sign dileptons plus jets with very little Standard 
Model (SM) background. With a NLSP gluino, however, the chargino as well as the leptons 
are absent in the gluino cascade decay. The conventional search strategy does not work here 
and the NLSP gluino can evade the current bounds from direct searches at the Tevatron. 
The NLSP gluino can be relatively light and its production rates can be large both at the 
Tevatron and especially at the LHC. 

On the other hand, as a consequence of t — b — r Yukawa unification, the third generation 
squarks, stops or sbottoms, are usually much lighter than those of the first two generations. 
In the large tan/3 limit, sbottoms are often the lightest squarks. As a consequence, gluino 
decays may lead to top-rich or bottom-rich events. In the co-annihilation region with AM = 
Mg — M^o <C 2rrit, there is no phase-space for on-shell top quarks. The gluino decay into 
6-jets, g — > bbxi, then becomes dominant and this enable us to search for NLSP gluino via 
multi-6 jets, namely 

pp,pp — > gg — > bbbb +0t- (1-2) 



Multi-6 events have been widely proposed for light scalar searches, for instance in NMSSM 
with the Higgs decaying into multi-6 jets via light scalars {h — > aa — > bbbb) fl~5[| . The multi-6 
jets plus significant missing transverse energy may only appear in Wh, Zh associated produc- 
tions whose production rates are much smaller than gluino pair production. The reconstruc- 
tion is also very different. 

The paper is organized as follows. In the next section (II), we discuss NLSP gluino 
decay and the parameter space where multi-6 jet production is significant. We very briefly 
comment models in which NLSP gluino can be realized and discuss its implications. In Section 
III and IV, we study the collider phenomenology of this scenario, with event selection and 
identification both at Tevatron and the LHC. We summarize our findings in Section V. 
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2. NLSP Gluino: Decays and Benchmark Models 



As mentioned earlier, the gluino-bino co-annihilation scenario requires the gluino to be NLSP 
in the sparticle spectrum, and to be nearly degenerate in mass with the bino LSP. The mass 
difference between the two should be < 20 %M^o . In the framework of minimal supergravity, 
this clearly requires non- universal gaugino masses at Mqijt- The leading motivations for 
grand unification theories (GUTs such as SO(10)) are the explanation of tiny neutrino masses 
and charge quantization. However, a partial unified model such as SU(A)a x SU(2) i x SU (2)# 
(4-2-2) also provides solutions to both problems. Non-universal asymptotic gaugino 
masses are naturally accomodated in the supersymmetric 4-2-2 model and have recently been 
investigated in Ref. pT| , |i~2^ ]. Other examples with non- universal gaugino masses include a 
supersymmetric SU{h) x 5J7(3)H ype rcoior proposed to explain doublet-triplet splitting problem 
in 517(5) GUT 0, [L8|, and GUT models with non-singlet F-term vevs @. 

Since we focus on a spectrum with NLSP gluino, on-shell charginos are kinematically 
forbidden in gluino decay. The color octet gluino therefore can decay only into colored SM 
particles such as the gluon octet or a qq pair, plus the color singlet LSP x?> 

g^qqXi,9Xi- (2-1) 

The three-body decay g — > qqx± is through off-shell squark interchange, while the two-body 
decay gxi can be realized by the triangle loop involving squarks. The partial widths of these 



two decay channels are given by [2C, 21] 
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Here Nib, Niyy and Nih u are respectively the bino B, wino W" and Higgsino H u components 
of the LSP neutralino Xi- The three-body decays will be suppressed if the scalar masses are 
too large, or due to phase space if the mass difference between g and Xi is too small. The 
two-body mode may dominate in this case. If either the gluino two-body decay g — > gxi, or 
three-body decay due to the small mass difference into light jets dominates, the final state jets 
are typically as soft as those from parton showers. In this case the gluino decay very likely 
gets buried in the huge QCD background, and consequently the search becomes extremely 
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challenging |]T^, p^j . Our study will focus on a different region where the gluino three-body 
decays into multi-6 jets dominate. 

To illustrate the search strategy, we have selected two benchmark points from a previous 
study of supersymmetric 4-2-2 models |TT| , O ], In 4-2-2 models, the matter fields of each family 
belong in (4,2,1) and (4,1,2). The third family fermion masses, to a good approximation, 
arise from a Yukawa coupling to the bi-doublet (1, 2, 2). Thus, t — b — r Yukawa unification 
arises as a natural prediction flTT| , 12, 14] 



Y t = Y b = Y T = y Dirac . (2.5) 

In SO(10), t — b — t Yukawa unification typically predicts gluino to be the lightest colored 
sparticle while in 4-2-2 models one realizes the NLSP gluino scenario through gluino-bino 
co-annihilation. A large bottom Yukawa also naturally drives sbottom to be the lightest 
squark. With AM ~ 50 GeV and O(TeV) squarks, g — >• bbxi decay often dominates. Figure |l] 
shows the dependence of the gluino decay BR in the AM — plane for the 4-2-2 model 
with /x < 0. We require consistency of the model with various phenomenological constraints 
such as BR(B S jiV), BR(6 s 7 ), BR(B tu), A(g - 2)„, WMAP relic density (in 
the 5(7 range), and all the sparticle mass bounds [11, 12]. In Fig. [l] the points which satisfy 



t — b — t Yukawa unification are a subset of the displayed points and mostly lie in the dense 
region around 1 TeV sbottom mass, where the mass difference AM is between 40 and 60 GeV. 
For this region the branching fraction of gluino three-body decay BR(<? — > bbxi) dominates 
both the two-body one, BR(g — > gxi), as well as the three-body decay into light quarks 
BR(g — > qqxi)- For a sufficiently large sbottom mass the two-body decay is dominant, as can 
be seen from Fig. |l]. Note that we show only those scenarios in Fig. |l| for which gluino is the 
NLSP. We have picked two benchmark models, one with ^ > (Model A), listed as point 1 
in Ref. [|llj; the second model has /i < (Model B), selected from the large number of points 
in Fig. [l]. The relevant observables are listed in Table |]. Both models can evade the direct 
search Tevatron gluino bounds due to the dominant decay mode g — > bt>Xi- 



3. NLSP gluino Search at the Tevatron 

In both benchmark models, the gluino masses are of order 200-300 GeV, and their pair 
production rates at the Tevatron are around the pico-barn level. Therefore, a search for a 
relatively light NLSP gluino appears quite promising. In this section, we illustrate how one 
could identify NLSP gluinos at the Tevatron for the two benchmark models above. 

As shown in Table |] for both models we focus on NLSP gluino with the dominant three- 
body decay g — > bb~Xi- Therefore we wish to identify signal events of gluino pair production 
with 4 6-jets plus missing transverse energy ^Ej-, 

pp -)• gg -)• bbbb+0 T . (3.1) 

Due to the relatively long lifetime of the S-mesons, their decays on average take place O(mm) 
distance away from the primary interacting vertex. With the vertex detector, tagging jet 
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Figure 1: Mass difference M g — M^o versus Mg j for 4-2-2 model with [i < 0. The points shown 
satisfy all the experimental constraints described in the text. The red circular points depict scenarios 
for which the branching fraction of two-body decay is dominant, i.e. BR(g — > gx\) > 40%. For the 
empty square points in black the three-body decay is dominant, BR(g — > bbxi) > 40%. Gluino is the 
NLSP for all the points shown in the figure. 



with decaying 5-mesons will significantly reduce the QCD jets background. The 6-tagging 
efficiency at the Tevatron is taken to be = 50% pi| ]. The 6-jet production in the SM is 
either due to gluon splitting or from top quark decay originating from top pair and/or single 
production (tt pair provides bb in the final states, while single top production jg — > btj' also 
provides bb), and so the 6-jets always arise as pairs. Hence we only need to tag three 6-jets so 
that we do not have to pay the additional 50% loss in the fourth 6-tagging. After multiplying 
by the 6-tagging efficiency, both the signal and the background events are reduced by 

4 = 12.5% . (3.2) 

By requiring 3 6-tagged jets, the SM production 46 + X becomes the leading irreducible 
background. Also, there exist reducible backgrounds due to other jets being mis-identified 
as 6-events. About 15% of D-mesons in the jets can be mis-identified as S's, and the mis-fr 



tagging rate for light jets is 0.4% [24]. Since the c-jet production in the SM is very similar to 
the 6-jet production, the production rate is at comparable level. With the 15% faking rate, 
we can safely neglect the c-jet production in our study. However, the production rate of light 
jets is several orders of magnitude higher than that of pure 6-jet production [25] and cannot 
therefore be neglected even with 0.4% mis-6-tagging rate. Therefore, we include jjbb + X in 
our background analyses. 
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Model A (jti > 0) 


Model B (/x < 0) 


m (GeV) 


14110 


1513 


Mi (GeV) 


499.54 


-479.49 


M 2 (GeV) 


832.03 


-845.5 


M 3 (GeV) 


0.7945 


69.53 


tan p 


50.82 


47.7 


A) 


-34551.2 


-1668.84 


m Hu (GeV) 


6092.74 


492.41 


m# d (GeV) 


14194.5 


1071.75 


M § (GeV) 


329 


261 


M^o (GeV) 


284 


207 


M Sl (GeV) 


5294 


950 


BR(g -> 66x?) 


76.3% 


50.8% 



Table 1: Model parameters at GUT scale (above double line) and low scale (below double line) for 
two benchmark models. Note that the bino component of Xi > 99.9%. 



Besides the bb, the dark matter particles Xi also appear in the gluino cascade decay. The 
missing transverse energy $t is another characteristic feature of the signal. The irreducible 
SM background for is from Z production, with the branching fraction of Z invisible decay 
(Z — > uv) as 20%. However, due to the uncertainty of mis-measurement in jet energy or 
momentum, the events without Z can also induce 0t- The third source is due to leptonic 
decays of bosons, especially — > t^v t — > i~^vi>y T v T where the leptons from r three- 
body decays are below the visible lepton cut {jf T > 10 GeV). The SM backgrounds that we 
consider in the study are then 

bbbbjjbb, 

bbbbZJjbbZ with BR(Z vv) = 20%, 

ti — > bbjjr^Vr with leptonic decay . (3.3) 

To simulate the detector effects, we smear the hadronic jet energy by a Gaussian distribution 
whose width is parameterized as [24] 

-FT 2 = , --, © 5%. 3.4) 

The following basic kinematical cuts on the transverse momentum (pt), the pseudo-rapidity 
(77), and the separation in the azimuthal angle-pseudo rapidity plane (A.R = \J (A(/>) 2 + (Ar]) 2 ) 
between two jets have been employed for jet-selection pi|1 

p> T > 15 GeV, 1^1 < 1.0 , ARjj > 0.4. (3.5) 
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Figure 2: (a) Minimal and maximal pr distribution of 4 6-jets in the signal pp — » gg — > bbbbxiXi a t 
Tevatron. (b) Missing transverse energy ( / E-r) distribution in signal events as well as its background 
AbZ and jjbbZ at Tevatron, taking the branching fraction of Z invisible decay (Z — > vv) as 20% and 
mis-6-tagging rate of light jet as 0.4%. 



Because of the relatively small mass difference AM = M g — M^o ~ 50 GeV, the 6-jets 
in final states could be rather soft. The minimal and maximal pt distributions of 6-jets are 
plotted in Fig. §(a). The plots that we show in this paper are all for benchmark model B, 
with the relevant features for Model A expected to be very similar. One can see that the 6-jet 
with minimal pt would be rejected by the basic pt cut in Eq. ( |3.5[) . To retain as many signal 
events as possible, we apply the pt cut for three 6-jets except the softest one. The softest jet 
will be viewed as unreconstructed calorimeter energy in the detector. It is consistent with 
the requirement of 3 tagged 6-jets above. The missing transverse energy J^t is reconstructed 
according to the smeared observed particles, namely 6-jets in our case. We show the $t 
distribution of the signal and backgrounds including basic cuts for rj and AR in Eq. ( |3.5| ) and 
6-jet and mis-6 jet tagging efficiency in Fig. |||(b). We find that by requiring a significant $t 
cut in the final states, namely 

Xt > 30 GeV, (3.6) 

the SM 4b and jjbb backgrounds at Tevatron can be completely eliminated. Since the signal 
events do not contain any lepton in the final states, we veto any event with visible leptons of 

p e T > 10 GeV . (3.7) 

The leading reducible background of soft leptons type is due to semi-leptonic tt events with 
one decaying into r ± , and further decaying into soft leptons. With one of the light 
jets faking a 6-jet, we find that the contribution due to tt is below 0.01 fb. The leading 
irreducible background after JSt cut is then jjbbZ with invisible Z decay. 
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c(fb)@ Tevatron 


Model A 


Model B 


46 


4bZ 


2j2bZ 


s/Vb 


basic cuts and 3b tagging 


2.3 


4.8 


2.7 x 10 3 


0.02 


1 




> 30 GeV 


1.4 


3.3 




0.019 


0.95 


4.5(A)/11(B) 



Table 2: Production cross section for pp — > gg — > bbbbxlxl m models A and B and for backgrounds 
46, AbZ, 2j2bZ, after basic cuts and missing energy cut at the Tevatron with a luminosity of 10 fb _1 . 



In Table |, we summarize the signal and background cross sections at the Tevatron for 
the two benchmark models, after implementing the basic cuts and cut. For our numerical 
analyses, we use the CTEQ6L1 parton distribution function |j26|| . The SM backgrounds are 
simulated by the automatic package Madgraph/Madevent |25[| . The signal significance is 
obtained in terms of Gaussian statistics, given by the ratio S/\fB of signal and background 



events with a luminosity of 10 fb . For benchmark model A with M g ~ 300 GeV, the 
statistical significance is close to 5 a, while for benchmark model B with a smaller gluino 
mass, as expected it is 2.5 times larger. 

In order to further identify the signal events and extract spectrum information, we pro- 
pose to reconstruct the events through the invariant mass distribution of di-6 jets M bb and 
Mt2 method. There are 4 6-jets in the final states so we will have three combinations of 6-jet 
pair in event reconstruction. Usually the two hardest 6-jets come from different gluinos and 
the 6-jets pair coming from the same g has smaller separation angles. Therefore, following the 
reconstruction method in Ref. ^3|, we select the two hardest 6-jets 61,62, and let 63 denote 
the third jet that minimizes A</>(62, 63), such that the pair 62, 63 come from the same g. The 
other 6-jets pair from the second g consists of two 6-jets 61 and 64. For gluino three-body 
decay, the dijet invariant mass M bb must be less than the mass difference between the gluino 
and the LSP masses [^7]] 

M bb < AM. (3.8) 

To illustrate this, we display the distribution of M bb = maxfMf,^, M b2b . a ] for benchmark model 
B in Fig. |3|(a); it clearly shows the edge of M bb around AM. It is important to note that 
the final states contain two invisible massive particles emanating separately from two parent 
particles, and therefore it is usually hard to reconstruct the signal events. However we still 
have sufficient information to fully reconstruct the signal events. Since the signal dominates 
according to our analyses above, the gluino mass can be estimated from the production rate. 
Once Mg and the mass difference AM are known, the LSP mass M^o can be easily obtained. 
By substituting the LSP mass M^o into the reconstruction, we can plot the variable Mt2 



which is defined as [28, [2 



M| 2 (g) = min {max[M^ (1) , Mj (2) ]}, (3.9) 

pf'+pf^pr 

where the transverse masses are constructed for two gluino decay chains in terms of 
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Tevatron Tevatron 




20 40 60 80 100 100 200 300 400 500 

M bb (GeV) M T2 (GeV) 

Figure 3: (a) The distribution of di-6 jet invariant mass M h ^. (b) Mti reconstruction with 200 GeV 
trial LSP mass. Both distributions are for Model B. 

the relevant transverse invariant mass and transverse momentum of bb system and as function 
of the trial LSP mass. Figure ^(b) shows the Mt2 reconstruction of the gluinos in Model B 
with trial LSP mass as 200 GeV. The consistency with our assumptions confirms the LSP 
mass M^o is the correct one and more importantly, it can also be used as a check for the 
gluino-bino co-annihilation condition as in Eq. (|l.l| ), 

M^o ~ Mg/(1 + 20%) ~ 200 GeV . (3.10) 
4. NLSP gluino Search at the LHC 

In this section we discuss the NLSP gluino search at the early LHC run of 7 TeV cm. energy. 
We employ a search strategy very similar to that for the Tevatron. In the LHC environment, 
the mis-6 tagging rate for light jets is about 1/30 for low pr range (15-50 GeV) |3(J. However, 
since the gluino pair production at the LHC is dominated by gg — > gg, the total cross section 
is about 30 times larger than at the Tevatron. We can then require to tag all four 6-jets and 
increase the missing transverse energy $t selection cut: 

• 4 b-tagged jets with pip > 15 GeV, |%| < 2.0, e& = 50% ; 

• ARjj > 0.4 ; 

• veto any event with lepton pj, > 10 GeV ; 

• T > 40 GeV . 



- 9 - 



cj(fb)@ 7 TeV LHC 


Model A 


Model B 


4b 


4bZ 


2j2bZ 


basic cuts and 4b tagging 


143 


271 


157 x 10 3 


0.55 


4.2 


> 40 GeV 


59 


140 




0.4 


3.3 



Table 3: Production cross section for signal pp —> gg — > bbbbxiXi an d backgrounds 4b,4bZ,2j2bZ 
after basic cuts and missing energy cut at 7 TeV LHC for the two benchmark models. 



The smearing parameterization is given as [3C] 

- 1 = : 3%. (4.1) 

Ej T^TGeV 

We show the px distribution of 6-jet in the signal final states and $t distribution for signal 
and backgrounds in Fig. ^a) and (b) respectively. One can see that the pr distribution of 
6-jets in the final states at the LHC is similar to that at the Tevatron. The missing transverse 
energy is smaller in comparison with the SM backgrounds. We summarize in Table | the 
results of signal and background studies at the LHC for the two benchmark models. After all 
the selection cuts, the signal events far exceed the SM background. By the first shut-down in 
winter 2010, the LHC should accumulate about 50 pb _1 of data. After the cuts, we expect 
negligible background events and 3 signal events for Model A, while Model B predicts about 
7 events. We expect our study to yield important clues about the underlying NLSP gluino 
scenario during the early stage of LHC operation. 

At the LHC, another channel may also become interesting. With a light gluino exchange 
in the i-channel, the first generation squarks u, d can be produced together with gluinos at 
significant rates via valence quark-gluon scattering. For Model B with Mq ~ ~ 1.5 TeV, 
the total production cross section for 

PP -> gu, gd -)• gu, gd (4.2) 

is about 120 fb. With a NLSP gluino, q —> gq decay always dominates and the gluinos will 
be highly boosted. The signal from heavy squarks decay will consequently be two extremely 
hard jets with one of them being the collimated gluino. It is then difficult to tag the 6-jet 
in decay products of boosted gluino. However, the 6-jets from the other gluino decay can 
still be tagged. The heavy squark resonance also provides a nice handle to suppress the SM 
backgrounds. In addition, with the heavy resonance, this channel may enable the search for 
g — > gxi decay. In either case, the search will require a more careful analyses of the jet 
substructure and we postpone this for future study. 



5. Conclusion 

We have explored the collider phenomenology of gluino-bino co-annihilation scenarios for 
both the Tevatron and 7 TeV LHC. The NLSP gluino is only slightly more massive (~ 
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Figure 4: (a) Minimal and maximal pr distribution of 4 6-jets in pp — > gg — > bbbb\iXi & t 7 TeV LHC. 
(b) Missing transverse energy (Xt) distribution in signal events and in backgrounds 4bZ and jjbbZ 
at 7 TeV LHC, with the branching fraction of Z invisible decay (Z —> vv) as 20% and mis-6-tagging 
rate of light jet as 1/30. 



50 GeV) than the bino LSP, so that the conventional gluino searches \\ X\ +jets are not 
applicable. We propose to search for gluino pairs through multi-b jets final states. By using 
two benchmark points from a supersymmetric 4-2-2 model, in which NLSP gluino arises 
naturally, we explicitly show how the search strategy works at hadron colliders. It is shown 
that with 10 fb _1 integrated luminosity, one can reach over 4.5 a at Tevatron. By the end 
of the first LHC run at 7 TeV with 50 pb^ 1 of accumulated data, the predicted signal events 
for both benchmark models are O(10), with negligible SM background events. 
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